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The polaron features for long-range electron-phonon interaction are in- 
vestigated by extending a variational approach previously proposed for the 
study of systems with local coupling. The ground-state spectral weight, the 
average kinetic energy, the mean number of phonons, and the electron-lattice 
correlation function are discussed for a wide range of model parameters focus- 
ing on the adiabatic regime and comparing the results with the short-range 
case (Holstein model). A strong mixing of electronic and phononic degrees of 
freedom for small values of the electron-phonon coupling constant is found in 
the adiabatic case due to the long-range interaction. Finally a polaron "phase 
diagram" is proposed. 
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I. INTRODUCTION 



In the last years many experimental results have pointed out the presence of strong 
electron- phonon {el —ph) coupling and polaronic effects in several compounds, such as high- 
temperature cuprate superconductors and colossal magnetoresistance manganites.^'^ This 
large amount of experimental data has renewed the interest in studying simplified el — ph 
coupled systems of the Holstein^ or Prohlich'^ type and in proposing more realistic interaction 
models. ^'^ 

The Holstein molecular crystal model is the prototype for short-range {SR) el — ph 
interaction since it takes into account local coupling of a tight-binding electron to optical 
local phonon modes. Till now an exact solution for this model has not been found and 
perturbative expansions are not able to describe the most interesting region characterized 
by intermediate el — ph couplings and electron and phonon energy scales not well separated. 
In this regime, as shown by several numerical studies^^^^ and variational approaches, ^^^^^ 
the system undergoes a crossover from a weakly dressed electron to a massive localized 
polaronic quasiparticle, the small Holstein polaron (SHP), with increasing the strength 
of interaction. All the ground state properties of the Holstein moldel can be described 
with great accuracy by a variational approach^^'^^ based on a linear superposition of Bloch 
states that describe weak and strong coupling polaron wave functions. Indeed this method 
provides an immediate physical interpretation of the intermediate regime and is in excellent 
agreement with numerical results. 

Recently a quite general el — ph lattice Hamiltonian with a "density-displacement" type 
interaction has been introduced in order to understand the role of long-range {LR) coupling 
on the polaron formation. ^'^^ The model for a single particle is described by the Hamiltonian 
H 

<i,j> i ^ i,j 

The units are such that h = 1. In Eq.(l) cj (q) denotes the electron creation (annihilation) 
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operator at site i, whose position vector is indicated by Ri, and the symbol <> denotes 
nearest neighbours hnked through the transfer integral t. The operator aj (oj) represents 
the creation (annihilation) operator for phonon on the site i, ujq is the frequency of the 
optical local phonon modes, a controls the strength of el — ph coupling, and f{\Ri — Rj\) is 
the interacting force between an electron on the site i and an ion displacement on the site j. 

The Hamiltonian ( 1 ) reduces to the Holstein model if / ( | — | ) = 5g, ^, , while in gen- 
eral it contains LR interaction. In particular when one attempts to mimic the nonscreened 
coupling between doped holes and apical oxygen in some cuprates,^ the interaction force is 
given by 

f{\R,-R,\)={\R,-R,\'' + iy\ (2) 

if the distance \Ri — Rj\ is measured in units of lattice constant. Considering the general 
el — ph matrix element 

M,= ^Y.fi\Rm\)e''-'^-. (3) 

with L number of lattice sites, we can define the polaronic shift E^, 

Ml 

^.-E-f, (4) 

and the coupling constant A = Ep/ zt, with z lattice coordination number, that represents a 
natural measure of the strength of the el — ph interaction in both SR and LR case. Clearly 
for LR interaction forces the matrix element is peaked around g = 0. Since it has 
been claimed that the enhancement of the forward direction in the el — ph scattering could 
play a role in explaining several anomalous properties of cuprates as the linear temperature 
behaviour of the resistivity and the (i-wave symmetry of the superconducting gap,^^'^^ the 
study of lattice polaron features for LR interactions is important in order to clarify the role 
of the el — ph coupling in complex systems. 

When the interaction force is given by Eq. (2), the model has been investigated applying a 
path-integral Monte-Carlo [PIMC) algorithm^'^^ efficient in the thermodynamic limit. The 



first investigations have been mainly limited to the determination of the polaron effective 
mass pointing out that, due to the LR coupling, the polaron is much lighter than the 
SHP with the same binding energy in the strong coupling regime. Furthermore it has 
been found that this effect due to the weaker band renormalization becomes smaller in 
the antiadiabatic regime. Then the quasi-particle properties have been studied by an exact 
Lanczos diagonalization method^^ on finite one-dimensional lattices (up to 10 sites) making a 
close comparison with the corresponding properties of HP. As a result of the LR interaction, 
the lattice deformation induced by the electron is spread over many lattice sites in the 
strong coupling region giving rise to the formation of a large polaron {LP) as in the weak 
coupling regime. All numerical and analytical results have been mainly obtained in the 
antiadiabatic and non-adiabatic regime. Only recently the behavior of the effective mass 
of a two-site system^^ in the adiabatic regime has been studied within the nearest-neighbor 
approximation for the el — ph interaction confirming that the LP is lighter than the SHP 
at strong coupling. 

In this paper we pursue the study of the ground state of the model with the interaction 
force given by Eq. (2) in the thermodynamic limit. We employ a variational approach 
previously proposed for the study of systems with el — ph local coupling^^'^^ and based on 
a linear superposition of Bloch states that describe weak and strong coupling polaron wave 
functions. Although the method is valid for any spatial dimension, we limit our study to the 
one-dimensional case. It has been found that the variational approach provides an estimate 
of the ground state energy in good agreement with PIMC results. The evolution of the 
ground-state spectral weight, the average kinetic energy, the mean number of phonons, and 
the electron- lattice correlation function with respect to the adiabaticity ratio uj^/t and the 
el — ph coupling constant is discussed focusing on the adiabatic regime. Indeed, in the 
adiabatic case, there is a range of values of the el — ph coupling where the ground state is 
well described by a particle with a weakly renormalized mass but a spectral weight much 
smaller than unity. Furthermore, with increasing the strength of interaction in the same 
regime, the renormalized mass gradually increases, while the average kinetic energy is not 
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strongly reduced. Finally regions of the model parameters are distinguished according to 
the values assumed by the spectral weight. The resulting "phase diagram" ^° shows strong 
mixing of electronic and phononic degrees of freedom for small values of the el — ph coupling 
constant in the adiabatic case. 



II. VARIATIONAL WAVE FUNCTION 

The variational approach is summarized following the lines of previous works. ^^'^^ 
We consider as trial wave functions translational invariant Bloch states obtained by 
taking a superposition of localized states centered on different lattice sites: 

\^f>=^i:e^'-^\i^l\Rn)>, (5) 



where 



(6) 



In Eq. (5) the apex i — w,s indicates the weak and strong coupling polaron wave function, 
respectively, |0 > denotes the electron and phonon vacuum state, and (f)^^\Rm) are variational 
parameters defining the spatial broadening of the electronic wave function. The phonon 
distribution functions h^^ik) are chosen in order to reproduce the description of polaron 
features in the two asymptotic limits. Therefore the weak coupling phonon distribution 
function h^y\k) is assumed as 

h^y\k)^ ^ ^, (7) 

where Eh{k) is the free electron band energy, while the strong coupling phonon distribution 
function h^^\k) as 

h'^\k)^^Y.\MRrn)\'e^'-^. (8) 

A careful inspection of weak and strong coupling wave functions shows that in the in- 
termediate regime they are not orthogonal and the off-diagonal matrix elements of the 



Hamiltonian are not zero. Hence the ground state energy is determined by considering as 
trial state a linear superposition of the weak and strong coupling wave functions: 



where 



\n / . . . ^ 14 / ... (10) 



and yS-g 



% = i^^L%^l^ (11) 



—(f) 

T 

is the overlap factor of the two wave functions > and iV'g >. In Eq.(9) Aj^ and 

are two additional variational parameters which provide the relative weight of the weak and 
strong coupling solutions of the system for any particular value of k. 

We perform the minimization procedure with respect to the parameters (j)^^\Rm), 
(t>^^\Rm), assuming 



S S' 5" 



where the quantities a^^\ (3^\ ^'^\ and rj^^^ denote variational parameters, and the symbols 
6, 5' , 5" indicate, respectively, the nearest, the next-nearest neighbors and so on. This choice 
takes into account the broadening of the electron wave functions up to third neighbors and 
provides an accurate description of the polaron features for any value of the parameters of the 
Hamiltonian. The ground state energies obtained with this choice are shghtly higher than 
PIMC mean energies, being the difference less than 0.5% in the worst case of intermediate 
regime. We note that these wave functions can be improved adding further terms in Eq. 
(12), so it is possible to obtain better and better estimates of the energy. 



III. RESULTS 

In this paper we study the properties of the ground-state in the one-dimensional case. 
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In Fig. 1(a) we report the polaron ground state energy for different values of the adia- 
baticity ratio as a function of the el — ph constant couphng a. We have checked that our 
variational proposal recovers the asymptotic perturbative results and improves significantly 
these asymptotic estimates in the intermediate region. Moreover, our data for the ground- 
state energy in the intermediate region are in very good agreement with the results of the 
PIMC approach^ shown as diamonds in Fig. 1(a). The consistency of the results with a 
numerically more sophisticated approach indicates that the true wave function is very close 
to a superposition of weak and strong coupling states. 

Another property of interest is the ground state spectral weight Z 

Z = Zk=o = \ < A=o\cUo\0 > \^ (13) 

that gives the fraction of the bare electron state in the polaronic trial wave function. It 
measures how much the quasiparticle is different from the free electron {Z — 1), and a small 
value of Z indicates a strong mixing of electronic and phononic degrees of freedom. As 
plotted in Fig. 1(b), the increase of the el — ph coupling strength induces a decrease of the 
spectral weight that is smooth also in the adiabatic regime. The reduction of Z is closely 
related to the decrease of the Drude weight obtained by exact diagonalizations^^ pointing out 
a gradual suppression of coherent motion. We note that the behavior of Z is different from 
that of the local Holstein model. In fact for the latter Z results to be very close to the ratio 
m/m*, with m and m* bare electron and effective polaron mass, respectively,^^ while for LR 
couplings Z < m/m* in the intermediate to strong coupling adiabatic regime. This relation 
is confirmed by the results shown in Fig.l (b), where the dash-dotted line and the squares 
on a similar line indicate the spectral weight Z and the ratio m/m* obtained within the 
variational approach, respectively, as a function of the coupling constant a at Wo/t = 0.25.^^ 
Actually there is a large region of the parameters in the adiabatic regime where the ground 
state is well described by a particle with a weakly renormalized mass but a spectral weight 
Z much smaller than unity. While the electron drags the phonon cloud coherently through 
the lattice, with increasing the el —ph coupling in the adiabatic case, a band collapse occurs 
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in the SR case, while the particle undergoes a weaker band renormalization in the case of 
LR interactions. Therefore in the LR case the polaron results lighter than the SHP in the 
intermediate to strong coupling adiabatic regime. 

Insight about the electron state is obtained by calculating its kinetic energy K in units 
of the bare one. Since the average kinetic energy gives the total weight of the optical 
conductivity, K includes both coherent and incoherent transport processes.^® As reported 
in Fig. 1(b) and Fig. 1(c), in the strong coupling adiabatic region before the electron is 
selftrapped {K <S 1), the average kinetic energy is weakly renormalized, the ratio m/m* is 
reduced and the spectral weight is nearly zero. Furthermore, in Figs. 1(c) and 1(d), within 
the adiabatic regime, the average kinetic energy and the mean number of phonons do not 
show any sharp change by increasing the el — ph coupling. 

Another quantity associated to the polaron formation is the correlation function S{Ri) 



or equivalently the normalized correlation function x{Ri) = S{Ri)/N, with N = J2iS{Ri). 
In Fig. 2(a) we report the correlation function S{Ri) at (Jo/t — 1 for several values of 
the el — ph interaction. The lattice deformation is spread over many lattice sites giving 
rise to the formation of LP also in the strong coupling regime where really the correlation 
function assumes the largest values. In the inset of Fig. 2(a) the normalized electron-lattice 
correlation function x shows consistency with the corresponding quantity calculated in a 
previous work.^^ While in the weak coupling regime the amplitude % is smaller than the 
quantum lattice fluctuations, increasing the strength of the interaction, it becomes stronger 
and the lattice deformation is able to generate an attractive potential that can trap the 
charge carrier. Clearly, even if the correlations between electron and lattice are large, the 
resulting polaron is delocalized over the lattice due to the translational invariance. Finally 
the variation of the lattice deformation as a function of cuo/t shown in Fig. 2(b) can be 
understood as a retardation effect. In fact, for small uJo/t, the less numerous phonons excited 
by the passage of the electron take a long time to relax, therefore the lattice deformation 





(14) 
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increases far away from the current position of the electron. 

In Fig. 3 we propose a "phase diagram" based on the values assumed by the spectral 
weight in analogy with the Holstein polaron.^"^ Analyzing the behavior of Z it is possible 
to distinguish three different regimes: (1) quasi-free-electron regime (0.9 < Z < 1) where 
the electron has a weakly renormalized mass and the motion is coherent; (2) crossover 
regime (0.1 < Z < 0.9) characterized by intermediate values of spectral weight and a mass 
not strongly enhanced; (3) strong coupling regime {Z < 0.1) where the spectral weight is 
neghgible and the mass is large but not enormous. We note that for LR interactions in the 
adiabatic case there is strong mixing of electronic and phononic degrees of freedom for values 
of the couphng constant A (sohd lines) smaller than those characteristic of local Holstein 
interaction (dashed lines). Furthermore in this case, entering the strong coupling regime, 
the charge carrier does not undergo any abrupt localization, on the contrary, as indicated 
also by the behavior of the average kinetic energy K, it is quite mobile. 

In order to study the effects of different el — ph interactions, we have evaluated the 
average kinetic energy for both LR and SR cases. As reported in Fig. 4(a), in excellent 
agreement with a previous study, for LR interactions K decreases very gradually with 
increasing A. Furthermore, if the regime of parameters where the spectral weight Z — 0.1 
is considered, as shown in Fig. 4(b), in the adiabatic case the average kinetic energy is 
larger for LR interactions (solid hne) with respect to local Holstein ones (dashed line) . The 
comparison emphasizes that due to LR interactions in the adiabatic regime K is slightly 
renormalized even if the coherent motion is small. 

IV. CONCLUSIONS 

In this paper we have extended a previous variational approach in order to study the 
polaronic ground-state features of a one dimensional el — ph model with long-range interac- 
tion. The trial function is based on a linear superposition of Bloch states that describe weak 
and strong coupling polaron wave functions and it provides an estimate of the ground state 
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energy in good agreement with numerical methods. The results relative to spectral weight, 
the average kinetic energy, the mean number of phonons, and the electron-lattice correlation 
function have been discussed mainly in the adiabatic regime. It has been possible to identify 
a range of intermediate values of the el — ph coupling constant in the adiabatic case where 
the system is well described by a particle characterized with a weakly renormalized mass 
but a small spectral weight. In the same regime, further increasing the el — ph coupling, the 
renormalized mass shows a smooth increase, while the average kinetic energy is not strongly 
reduced. Finally we have proposed a "phase diagram" according to the values assumed by 
the spectral weight. It is found that there is strong mixing between electronic and phononic 
degrees of freedom for small values of the el — ph coupling constant in the adiabatic case. 

The variational approach can be easily generalized to high dimensions, and it has been 
recently applied to study the three dimensional continuum Frohlich model giving a very 
good description of ground state features. In any case the results discussed in this paper 
are not limited to the one-dimensional case as confirmed by the behavior of some properties 
on the square lattice. ^'^^ 

FIGURE CAPTIONS 

Fig.l The ground state energy Eq in units of loq (a), the spectral weight Z (b), the average 
kinetic energy K (c) and the average phonon number N (d) as a function of the 
coupling constant a for different values of the adiabatic ratio: cuo/t — 2 (solid line), 
ujo/t — 1 (dashed line), uo/t — 0.5 (dotted line) and too/t — 0.25 (dash-dotted line). 
The diamonds in figure (a) indicate the PIMC data for the energy kindly provided 
by P. E. Kornilovitch uq/I — 1, and the squares on a dash-dotted hue in figure (b) 
denote the ratio m/m* obtained within the variational approach at uj^/t — 0.25. 

Fig.2 (a) The electron-lattice correlation function S{Ri) at Uq/I — 1 for different values 
of the coupling: A = 0.5 (circles), A = 1.25 (squares), A = 2.0 (diamonds), A = 2.75 
(triangles up), and A = 3.5 (triangles down). In the inset the normalized correlation 
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function xi^i) ^'^ ^o/t = 1 for A = 0.5 (circles) and A = 2.75 (squares). 

(b) The electron-lattice correlation function S{Ri) at a = 2 for different values of the 
adiabatic parameter: cuo/t — 2 (circles), cuo/t — 1 (squares), cuo/t — 0.5 (diamonds), 
and LOo/t — 0.25 (triangles up). 

Fig. 3 Polaron "phase diagram" for long-range (sohd line) and Holstein (dashed line) el—ph 
interaction. The transition lines correspond in weak coupling to model parameters such 
that the spectral weight Z = 0.9, in strong coupling such that Z = 0.1. 

Fig. 4 (a) The average kinetic energy as a function of the constant coupling X at ujo/t — 1 
for long-range (solid line) and local Holstein el — ph interaction (dashed line). 

(b) The average kinetic energy X as a function of the adiabaticity ratio cuo/t for long- 
range (solid line) and Holstein (dashed line) el — ph interaction in correspondence of 
model parameters such that the spectral weight Z — 0.1. 
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